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ABSTRACT 

We study the rms peculiar velocity of clusters as a convenient statistic to put con- 
straints on cosmological models. This statistic is easy to compute theoretically given a 
model for the power spectrum; we show that with some assumptions it can be directly 
related to the quantity Hq' 6 ^. We develop a likelihood method for estimating the rms 
peculiar velocity of clusters from line-of-sight velocity measurements and their associ- 
ated errors. We apply our method to two samples of cluster peculiar velocities; a new 
sample known as the SCI sample and a subsample of the Mark III catalog. Although 
these two samples initially give results which are inconsistent, we show that they can 
be put into good agreement by the removal, particularly from the Mark III sample, 
of a set of clusters with multiple mass concentrations, cluster membership ambigui- 
ties, or excessive obscuration, properties which may have introduced unaccounted for 
errors into measurements of their peculiar velocities. Once these clusters are removed 
from the samples, they both favor a relatively low value for the 1-D rms peculiar ve- 
locity of clusters a v = 265^75 6 km s _1 (at 90% confidence), leading to the constraint 
SIqVs = 0.441° consistent with cluster abundance studies but inconsistent with 
fi = 1 CDM normalized to COBE at the 99.7% confidence level. 
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1 INTRODUCTION 

Several recent papers have shown that the peculiar veloci- 
ties of clusters of galaxies can be used to put strong con- 
straints on cosmological parameters (Croft and Efstathiou 
1994; Bahcall, Gramann, and Cen 1994; Bahcall and Oh 
1996; Moscardini et al. 1996; Borgani et al. 1997; Shi 1997). 
Much of their utility in this regard derives from the fact 
that the motions of clusters probe the peculiar velocity field 
predominantly on scales where it is still linear, making the 
comparison to theoretical models relatively straightforward. 
In addition, peculiar velocities of clusters can in principle 
be measured to much higher accuracy than those of galax- 
ies, since many galaxies can be measured per cluster. Thus 
the motions of clusters should provide a valuable tool in the 
study of large scale structure. 

However, results obtained thus far using different cat- 
alogs of cluster peculiar velocities have been contradictory 
(see, e.g., Moscardini et al. 1996). In order to have confi- 
dence in the use of cluster peculiar velocities as a probe 
of large scale structure it is important to understand why 
this is so. In this Letter we study this issue, focusing on 
the estimation of the rms peculiar velocity of clusters. This 
statistic gives us a convenient way to compare different sam- 
ples of cluster peculiar velocities and can be compared di- 



rectly to predictions of cosmological models obtained from 
either iV-body simulations or linear theory. Indeed, as we 
show below, the rms peculiar velocity of clusters can be 
directly related to the combination of parameters fig 6(J 8, 
where <rg is the rms mass fluctuation on a scale of 8 /i -1 Mpc 
(#o = 100 h km s- 1 Mpc^ 1 ). 

We use a likelihood analysis to estimate the rms pecu- 
liar velocity of clusters from line-of-sight velocities and their 
estimated errors. This analysis allows us to assess the consis- 
tency of two samples in a model independent way. We apply 
our method to two well-studied catalogs of cluster peculiar 
velocities; a new sample based on measurements of 782 spi- 
ral galaxies in the fields of 24 clusters described in Giovanelli 
et al. (1997) (hereafter the SCI sample) and a subset of the 
cluster peculiar velocities contained in the Mark III catalog 
(Willick et al. 1995, 1996, 1997a) based on measurements of 
346 galaxies in the fields of 31 clusters described in Mould et 
al. (1991), Han (1992), and Mould et al. (1993) (hereafter the 
HM sample) . The HM sample was used to set the zero-point 
of the Tully-Fisher relation for the Mark III catalog. These 
samples cover the sky uniformly, have comparable depths, 
and are based on Tully-Fisher distance determinations. As 
we shall show below, these two samples yield results which 
are inconsistent at a high level. 
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We argue that the inconsistency between the two sam- 
ples is due to the inclusion of clusters, mostly in the HM 
sample, whose associated galaxy samples are not represen- 
tative of a single collapsed object or where measurement 
difficulties make error estimates unreliable. We show that 
when clusters suspected of being in this class are removed, 
the SCI and HM samples show excellent agreement, yielding 
a value of the rms peculiar velocity of clusters ~ 250 km s _1 . 
This favors values of f2o' 6 U8 ~ 0.5, in agreement with cluster 
abundance studies. Our results emphasize that the careful 
selection of the galaxy sample used in the analysis of a clus- 
ter is crucial for the accurate determination of its peculiar 
velocity. 



2 ANALYSIS 

Linear theory of the growth of perturbations suggests that 
the I-D peculiar velocity dispersion a 2 = (vf D ) of clusters 
is given by an expression of the form 



(1) 



i r7-2 1.2 r°° 

^ = \^-j P(k)W 2 (k)dk, 

where P(k) is the power spectrum and W 2 (k) is a window 
function that filters out power on scales smaller than a clus- 
ter size. Croft and Efstathiou (1994), Bahcall et al. (1994), 
and Colberg et al. (1997) have shown that Eqn. 1 gives good 
agreement with N-body simulations when a Gaussian win- 
dow function is used, W 2 {k) = exp(~k 2 R 2 ), although they 
disagree on the best value for the parameter R. We choose 
to adopt the value 72 = 1.5 h~ Mpc used by Colberg et al. 
(1997), which works well for models normalized to cluster 
abundances. Choosing the larger value of 3.9 /i _1 Mpc used 
by Borgani et al. 1997 only changes our results by ~ 15%. 

Assuming a scale- invariant spectrum of initial mass fluc- 
tuations, we follow Efstathiou, Bond, and White (1992) and 
parametrize the power spectrum as 

P(u) = Eh (9) 

y> {l + [afe + (6fc) 3 /2 + ( c fc)2]-}2/^' ^' 

where a = (6.4/r) ft^Mpc, b = (3.0/r) ft^Mpc, c = 
(1.7/r) /i _1 Mpc, and v = 1.13. Here the parameter F 
(= Qoh in models with flo + A = 1) determines the "shape" 
of the power spectrum. Although motivated by CDM mod- 
els, this parametrization is able to fit a a wide range of 
theoretical models, including those that include mixed dark 
matter. 

The normalization B can be expressed in terms of the 
rms linear mass fluctuation on scales of 8 h~ Mpc, as, which 
is defined by 



2 
^8 



P(k)Wg{k) d 3 k, 



... . ■ / ■ - (3) 
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where W s (k) = 3ji (kR s )/kR 8 with R 8 = 8 7i _1 Mpc and ji 
is a spherical Bessell function. Following Colberg et al. 1997, 
we can express a v as 

100 km s" 1 o.f 



V3 



ng-V 8 v7(f), 



where the function f(T) is defined as 
/ P(k)W 2 (k) dk 



f(T) = 4nh' 



J P(k)Wi(k) d 3 k' 



(4) 



(5) 



and is fit well by the equation f(T) = 12. 5r -1 - 08 +49.4. Thus 
we can place an upper bound Q, ( Q 6 a$ < a„(405 km s -1 ) -1 . 
Alternatively, we note that v7 varies only between 11.0 and 
8.7 over the range of values F = 0.2 — 0.5. If we make the 
reasonable assumption that F lies in this range (see, e.g., 
Peacock and Dodds 1994), then to a good approximation 
77 = 10 and n[]' 6 (j 8 = ct„(580 km s" 1 )" 1 . Colberg et al. 
(1997) take a different approach, noting that if one adopts 
the value of as obtained from cluster abundance studies then 
the quantity Q,Q A 'as varies only very weakly with Qo, so that 
a v is nearly a function of F alone. Here we choose to view 
the measurement of a v as an opportunity to independently 
verify the results of cluster abundance measurements. 

A problem that arises in calculating a v from cluster 
data is that even if actual peculiar velocities are gaussian 
distributed, measured velocities are not. This is due to the 
fact that each velocity measurement has its own distinct er- 
ror. For distant objects, these errors can often be comparable 
to the velocities themselves and must be properly accounted 
for. Previous studies have relied on Monte Carlo methods to 
generate ensembles of simulated data sets convolved with 
measurement errors, which are then compared to the actual 
data. Here, we introduce a likelihood analysis by which the 
rms peculiar velocity of clusters can be estimated directly 
from measured line-of-sight velocities and their errors. 

We assume that each measured line-of-sight velocity w» 
is drawn from a gaussian distribution with variance a 2 + a 2 , 
where at is the measurement error in the measurement of 
Vi. With this assumption, the likelihood of measuring a set 
of Vi given a particular value of a v and measurement errors 
ai is given by 



L(a v ) = Yl 



(27^/2(^+^)1 



/2 



exp 



2(a 2 



(6) 



Bayes' Theorem tells us that the posterior probability 
P(a v \vi,ai) = L(a v )P(a v ), where P(a v ) is the prior and 
contains any information that we might know about a v be- 
forehand. In this paper we shall make the simplest assump- 
tion that P(a v ) = constant, so that P(a v \vi,ai) oc L(a v ). 
For a given set of velocities and errors (vi,ai), Eqn. 6 can 
be used to calculate the maximum likelihood value of a v as 
well as confidence intervals. We shall choose confidence in- 
tervals for a v using the criterion that L(a v ) > C, where C 
is determined by the desired confidence level. 

Models of large scale structure can be accepted or re- 
jected by comparing their predicted value of the peculiar ve- 
locity dispersion, a prc d to the likelihood distribution L(a v ). 
In particular, a confidence level for the rejection of a model 
can be assessed by calculating the relative area under the 
likelihood function satisfying L(a v ) < L(a plc<1 ). 



3 RESULTS 

In Fig. 1 we show the likelihood function L(a v ) for the SCI 
and HM cluster samples, both normalized so that they have 
unit area. The relatively low value of a v favored by the SCI 
sample is consistent with the value of a v ~ 250 km s _1 
predicted by models normalized to reproduce the correct 
abundance of clusters (Colberg et al. 1997), and is in- 
consistent at the 99.8% confidence level with the value of 
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Figure 1. The normalized likelihood functions for a v obtained 
from the SCI and HM samples. 



Figure 2. Peculiar velocities and lcrerrors of the 12 objects com- 
mon to the SCI and HM samples. 



<T V = 489 km s _1 predicted by an Q,o = 1 CDM model nor- 
malized to microwave background fluctuations (Bahcall et 
al. 1994). (Note that these papers quote the 3-D rms ve- 
locity, which must be divided by \/3 for comparison to a v , 
the 1-D rms velocity.) The larger value favored by the HM 
sample is quite consistent with £7o = 1 CDM normalized to 
COBE fluctuations. 

We can quantify the disagreement between the two sam- 
ples using their normalized likelihood functions. The prob- 
ability of the two samples having values of a v which differ 
by A is given by P(A) oc J o °° L H m(A + a v )L S ci(<?v) da v . A 
measure of the probability that the two data sets are drawn 
from distributions with the same a v is given by the relative 
area under the curve satisfying P(A) < P(A — 0). Using 
this measure we find that the SCI and HM samples disagree 
on the value of a v at the 99.8% confidence level. 

Insight into the nature of the discrepancy between the 
two samples can be gained by examining the 12 objects that 
they have in common. In Fig. 2, we show the velocity and 
error estimates for these 12 clusters as given in the SCI (x- 
axis) and HM (y-axis) catalogs. The measurements are in 
good agreement except for the Cancer cluster, corresponding 
to the point in the lower middle of the plot. This cluster has 
been studied by Bothun et al. (1983), who concluded that it 
consisted of several discrete groups spread in velocity over 
2800 km s _1 . Furthermore, their results suggest that the 
groups are not gravitationally bound, but rather are separat- 
ing with the Hubble flow. Giovanelli et al. (1997) selected all 
of their Cancer cluster galaxies from the largest mass group, 
whereas the Cancer cluster galaxies in the HM sample are 
drawn from several of the groups. The good agreement of 
the reported velocities for the less complex clusters suggests 
that the inconsistency regarding the Cancer cluster is due to 
the difference in membership criteria. Indeed, it seems likely 



that the velocity from the HM sample does not reflect that 
of a single gravitationally bound object. 

In order to explore the possibility that the discrepancy 
between the SCI and the HM samples is due to effects that 
are not properly accounted for in the error estimation, we 
construct subsamples SCI, and HM, of the cluster samples 
from which we remove clusters that we consider suspect. 
These mostly consist of clusters which the observers noted as 
having ambiguities in their associated galaxy samples or par- 
ticularly poor data. We also remove clusters for which other 
observers have noted evidence for multiple mass concentra- 
tions or overlap with other clusters. From the SCI sample, 
we excise six paired clusters (A2197/A2199, Pavo/Pavo II, 
and A2634/A2666) which have possible cluster membership 
ambiguities. A2634 and Pavo II are also in the HM sample, 
and we remove them from this sample as well. The following 
additional clusters are removed from the HM sample: Can- 
cer (see above), Pisces (Sakai et al. 1994), MKW1S, and 
A779 (Mould et al. 1993), due to their having multiple mass 
concentrations; 3C 296, due to its distribution of galaxies on 
the sky not being well defined (Mould et al. 1993); 1559+19, 
due to it being sparsely sampled and composed of rather in- 
homogenious data (Mould et al. 1993); OCA 3627, due to 
difficulties arising from its low galactic latitude (Mould et 
al. 1991); N5419, due to poor sampling and possible sub- 
structure (Mould et al. 1991); and OCA 3560 (Mould et al. 
1991), ESO 508, and CEN 45 (Aaronson et al. 1989), due to 
possible contamination by foreground or background galax- 
ies. This leaves us with reduced samples of 18 SCI* clusters 
and 18 HM, clusters. 

In Fig. 3 we show normalized likelihood functions for a v 
obtained from the SCI, and HM* subsamples. While the re- 
sults for the SCI sample are essentially unchanged, we see a 
marked reduction in the maximum likelihood value of a v for 
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Figure 3. The normalized likelihood functions for a v obtained 
from the SCI* and HM, samples after the removal of suspect 
clusters. 



the HM, clusters. We have determined that the probability 
of selecting a subset of 18 clusters at random from the HM 
sample that results in a maximum likelihood value of a v as 
small or smaller than that of the HM, sample is only 0.002. 
This strongly suggests that the suspect clusters in the HM 
sample have systematically higher measured velocities rela- 
tive to their reported errors than the other clusters in the 
sample. Given that the excised clusters were chosen due to 
ambiguities in their galaxy samples, it seems likely that these 
ambiguities have resulted in contributions to their measured 
velocities from unaccounted-for errors. The fact that the SCI 
results are not changed significantly by the removal of the 
suspect clusters is probably due to the strict membership cri- 
teria employed by Giovanelli et al. (1997) in the selection of 
their cluster galaxy samples. The use of these strict criteria 
lowers the probability of incorrect membership assignment 
and helps to ensure that each galaxy sample describes a sin- 
gle collapsed object, thus reducing the possible errors in the 
velocities measured for the 6 suspect clusters. 

With the suspect clusters removed, the two subsamples 
are in good agreement, both favoring a relatively low value 
of a v . The SC* subsample, which has smaller errors than the 
HM, subsample, yields a v — 247lgQ 3 km s _1 (here and be- 
low errors are quoted at the 90% confidence level). Making 
the assumptions about the power spectrum outlined above, 
this gives f2o' 6 <Tg = 0.421°' 14, in good agreement with a pre- 
vious analysis of the SCI data by Borgani et al. (1997). We 
get only a slight improvement in the accuracy of our re- 
sult by combining the SC* and HM, samples; keeping the 
measurement with the smaller error for clusters present in 
both results in a combined sample of 28 clusters yielding 



consistent with the value of 489 km s predicted for stan- 
dard Q = 1 CDM normalized to COBE at the 99.7% con- 
fidence level. 



4 DISCUSSION 

We have shown that when suspect clusters are removed, 
the SCI and the HM samples show excellent agreement, 
both favoring relatively low values for the rms peculiar ve- 
locity of clusters; using a combined sample we find a v — 
2571 re 8 km s _1 . If one makes some reasonable assumptions 
about the shape of the power spectrum, this converts to 
fio 6 <7 8 = 0.44±g;^|. These results are in good agreement 
with the results of cluster abundance studies, but are incon- 
sistent with Q = 1 CDM normalized to COBE at the 99.7% 
confidence level. 

Our results suggest that peculiar velocities calculated 
for galaxy samples that are not representative of a single 
mass concentration are susceptible to large unaccounted- 
for errors. This is perhaps not too surprising, since veloc- 
ities are calculated assuming that all of the sample galaxies 
have a common distance. While this assumption is always 
violated to some extent due to the non-zero thickness of a 
cluster, the presence of multiple mass concentrations or con- 
tamination by foreground or background galaxies makes it 
much more likely that the galaxy sample will by distributed 
anisotropically about the cluster center, leading to a bias in 
the measured velocity. This problem is especially severe for 
spiral galaxy samples, since they tend to inhabit the outer 
regions of clusters. The fact that the central cores of clus- 
ters consist primarily of ellipticals makes it likely that as the 
fundamental-plane distance measure continues to improve it 
will become the preferred method for measuring the peculiar 
velocity of clusters. 

It is unclear what implications, if any, these results have 
for the Mark III catalog of galaxy peculiar velocities. Al- 
though the HM cluster data set was used in the calibration 
of this catalog, Willick et al. (1995) recognized that not all 
of the galaxy samples were well modeled as belonging to 
a single cluster and accounted for this in their calibration 
procedure. 

More cluster peculiar velocities are needed in order to 
improve the robustness of our results. In particular, a larger 
sample will allow us to study the dependence on cluster 
selection and the validity of our approximation that cluster 
peculiar velocities are gaussian distributed. However, our 
results underscore the importance of careful galaxy sample 
selection when measuring cluster peculiar velocities. 



a v = 257±^ 8 km s" 1 and fig- 



3 a s = 0.44±°;ijj. Using the 



likelihood function we find that the combined sample is in- 
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